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ABSTRACT 

Cold, dense clouds of gas have been proposed as baryonic candidates for the 
dark matter in Galactic haloes, and have also been invoked in the Galactic disc 
as an explanation for the excess faint sub-mm sources detected by SCUBA. 
Even if their dust-to-gas ratio is only a small percentage of that in conven- 
tional gas clouds, these dense systems would be opaque to visible radiation. 
This presents the possibility of detecting them by looking for occultations of 
background stars. We examine the possibility that the data sets of microlens- 
ing experiments searching for massive compact halo objects can also be used 
to search for occultation signatures by cold clouds. We compute the rate and 
timescale distribution of stellar transits by clouds in the Galactic disc and 
halo. We find that, for cloud parameters typically advocated by theoretical 
models, thousands of transit events should already exist within microlensing 
survey data sets. We examine the seasonal modulation in the rate caused by 
the Earth's orbital motion and find it provides an excellent probe of whether 
detected clouds are of disc or halo origin. 



1 INTRODUCTION 

Over the last decade several workers have suggested that the Galaxy may contain a large 
population of cold (T < 10 K), dense, low-mass clouds of molecular material. Several authors 
have argued that the mass contained in such clouds could be responsible for keeping galaxies' 
rotation curves flat beyond the edges of optical discs. Pfenniger fc Combes (199^ argued 



that the clouds would be an integral part of a fractal interstellar medium and be confined 
to a thin disc. By contrast, de Paolis et al. (1995 ), Gerhard fc Silk (1996 ), Draine (1998 ), 
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and [Walker fc Wardle (1998| ) considered the case of cold clouds that are distributed in an 
approximately spherical halo. De Paolis et al, Draine and Gerhard & Silk argued for clouds 
of approximately solar mass, while Walker & Wardle pointed out that ionized winds from 
clouds with masses around a Jupiter mass ( Mj) could explain the extreme scattering events 
that are frequently observed in the light curves of pulsars seen at low Galactic latitudes. 

A major issue with these proposals that a significant mass is locked up in very small 
clouds, is to understand why the clouds do not collapse to form stars or degenerate objects 
such as brown dwarfs. Sciama (200(j| ) has argued that cosmic-ray heating may be able to 
balance molecular cooling in a stable way for dust-free clouds, whilst [Lawrence (2001[ ) has 
presented similar arguments that cosmic-ray heating could balance cooling by dust emission 
in dusty clouds. Clouds that are an integral part of the ISM would be expected to be dusty, 
and radiate strongly at sub-mm wavelengths, p^awrence (2001| ) has considered the possibility 
that a significant fraction of the reddest SCUBA sources could lie in the Galaxy. He concludes 
that for dusty clouds, masses between ~ 10 Mj and ~ 10 Mq can be ruled out, and that 
a population that contributes significantly to the SCUBA counts will make a significant 
contribution to the mass of the Galaxy only if they have masses ~ Mj. These objects would 
have to be confined to the Galactic plane. 

The clouds would have temperatures T ~ 10 K, and if they were in virial equilibrium, 
they would be ~ 20 AU in diameter and totally opaque at optical wavelengths (and even 
have an optical depth of order unity at 1 mm). Their number density on the sky has to be 
comparable to that, ~ 1000 deg~^, of unidentified SCUBA sources, so they would cover a 
fraction ~ 10~^ of the sky. Consequently, at any given time, of order one star in a million 
should be occulted by one of these clouds. It should be feasible to detect such occultations 
in the databases that have been produced in searches for microlensing events (for a recent 
overview of microlensing results see [Kerins 200 1| , and references therein). 

Here we investigate the rate at which stars would be occulted by opaque clouds. We 
argue in Section ^ that halo and disc clouds are likely to have the characteristics required 
to produce observable transits. In Section ^ we derive general formulae for the rate and 
timescale distribution of transit events and in Section ^ we apply them to a specific Galactic 
model. Currently available microlensing databases refer to a handful of lines of sight. The 
largest body of data refers to lines of sight in the general direction of the Galactic centre. 
Other extensive data sets refer to sight lines towards the Magellanic Clouds. Consequently, 
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we concentrate on predictions for these particular directions. We also investigate means by 
which halo and disc cloud populations could be differentiated. 



2 CLOUD CHARACTERISTICS AND DETECTABILITY 

In order to occult background stars cold clouds must meet certain criteria. Firstly, the clouds 
must be sufficiently diffuse that they do not appreciably gravitationally lens the starlight. 
The case of gravitational lensing by gas clouds has been studied previously by [Henriksen fc 



Widrow (f995| ). A necessary condition to avoid the gravitational lensing regime is that the 



cloud radius R exceeds the Einstein radius, leading to the requirement that 

where M is the cloud mass and D = D^^Ds — Dc)/Ds, with (-Dg) the distance between 
the observer and cloud (source). 

The long-term stability of cold (T ~ 10 K) clouds is an open issue and a serious the- 
oretical problem for the scenario. If opaque clouds can be stabilized against collapse, their 
size would be of the same order as their virial radius (c.f. [Lawrence 200 1| ) 



R... = 7AV( f ](^) . (2) 
\10-^MqJ VIO KJ ^ ' 

R ~ -Rvir would also satisfy the lensing condition in equation (|l|) provided 

f T \^ ( D 



M»h6xl0-Me^^j ^^j. (3) 

This condition is satisfied for almost all plausible clouds since, regardless of their formation 
mechanism, hydrogenous clouds less massive than ~ 10^^ Mq evaporate away over a Hubble 
time (|de Rujula, Jetzer fc Masso 1992|) . 



Diffuse clouds will still refract, rather than occult, starlight if they are transparent 
( [Walker fc Wardle 1998| ; Praine 1998[) . How opaque might we expect the clouds to be? 



iVH = 3.3 X 10- m- (t^) • (4) 



The average hydrogen column density through a virialized cloud of is 

R / T 

Assuming typical dust properties, the visual extinction through a cloud is therefore (e.g. 
Binney fc IVlerriheld 1998[ ) 
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where / is the factor by which the cold cloud dust-to-gas ratio exceeds that of typical 
Galactic molecular clouds. Equation implies that if cold clouds have even remotely the 
same dust abundance as normal clouds then they easily fulfill the requirement that they be 
opaque to starlight. In fact, for the clouds not to produce a detectable occultation (which, 
conservatively, requires Ay < 1) they must have a dust-to-gas ratio at least four orders of 
magnitude lower than is measured for ordinary gas clouds, i.e. they must be essentially dust 
free. 

We should expect disc clouds to have a similar dust-to-gas ratio as the interstellar medium 
and therefore be completely opaque. However, equation (j^) also implies strong limits on halo 
cloud properties for them to escape detection by occultation. This can be more clearly seen 
by considering how halo clouds may have interacted with the interstellar medium over the 
lifetime Tgai of the Galaxy. Over this time one would expect a typical halo cloud to have 
passed through the Galactic disc times, where 



a 



ll60kms-ij 1^10 GyrJ I^SkpcJ ' 
with Tc and a the halo core radius and velocity dispersion. Neglecting accretion the mass 
swept up by these clouds on each passage through an exponential disc is 

R 



7 AU 

^ (O.O3M0PC-3) (200 pc) ' ^'^^ 
where p is the typical density of gas in the plane of the disc, H the disc scale height and e 
the efficiency with which the clouds retain swept up gas. Adopting the parameter values in 
equations (|^) and implies that the present mass fraction of swept-up disc gas in a typical 
halo cloud is NMs^cep/M ~ O.OOTe for a 10~^Mq cloud. This, together with equation (^, 
indicates that provided e > 0.008/^^ the cloud dust-to-gas ratio should be large enough to 
produce at least a one magnitude diminution in the light from background stars. 

Lastly, we need to know how easily the flux change from transit events can be detected. 
An important consideration is source resolvability. For Galactic microlensing experiments it 
is often the case that more than one star occupies the seeing disc, leading to a detection 
bias referred to as blending (|Alard 1997 ). The situation is more extreme for experiments 
looking towards M31 where there are many stars per detector pixel. Transit events could 
escape detection by these surveys if the occulted star lies so close to other stars that the 
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effect of its occultation is masked by light from otlier stars. Quantitatively, if a fraction of 
the flux in the seeing disc comes from the target star, the flux at minimum will be a fraction 
(1 - /*) + /*10-°-^^^ of the flux prior to occultation. Equating this fraction to iq-o-^^^t^ 
where AMt is the smallest magnitude change a survey can measure, we find that the event 
will be detected only if 



For completely opaque clouds (large Ay) the requirement is simply > AMp for AMp <S 1. 
In the context of gravitational microlensing searches, experiments to detect microlensing 
events in our Galaxy typically adopt the very conservative value AMp ~ 0.3 as a means of 
excluding variable stars. The M31 surveys typically use a much lower threshold AMp ~ 0.02. 
If similar thresholds are adopted for the detection of transit events then they should be 
detectable within our Galaxy provided blending effects are not too severe. It may even 
be possible to detect transit events towards M31 provided the occulted stars ordinarily 
contribute upwards of 2% of the flux in a seeing disc. 

3 TRANSIT RATE AND TIMESCALE DISTRIBUTION 

From here onwards we assume that the cold clouds are opaque and exceed their Einstein 
radius, so that they give rise to observable occultation events. We calculate the rate at which 
occultations of a given duration should be observed at any season of the year and along any 
hue of sight. Similar calculations of microlensing rates are described by Griest (1991) and 
Kiraga & Paczyhski (1994). 

Let there be n(v)d^v spherical clouds of radius R per unit volume and with velocities 
in the plane of the sky in d^v. Then in a sheet of thickness dD^ that extends perpendicular 
to the line of sight, the rate at which clouds pass a given background star with impact 
parameters in the range (6, b + db) is 




1 - 10-0-4^^ 



(8) 



d^T = n{v)d\dD^2v±db, 



(9) 



where v± — |v|. In an obvious notation we have d^v = v±dv±d9, so 



d^r = 2n{v)vldv± dOdD^db. 



(10) 



An impact with b < R causes an occultation of duration 



t = 2VR^-b^/v±, 



(11) 
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so we may obtain the distribution over durations by eliminating h in favour of t. One finds 

n{y)v\dv^ dO 



d^r 



(12) 



d^cdt ^(2R/vi_tY - 1 
This equation gives the occultation rate of a single source star. To get the rate that can be 
measured observationally, we have to average it over all detectable source stars (c.f. [Kiraga 
fc Paczyhski 199^ ). If the probability that a source star at distance has its velocity on 
the sky in d^Vg is p(vs)d^Vs, then the required average of equation ([12D is 
d^r /■ ,0 . . r^ij/t 



dDAt 



/f 



{2R/v±ty 



(13) 



v± is the magnitude of the cloud's transverse velocity with respect to its local point on the 
line of sight from the source to the Sun; consequently it is a function of Vg. The argument 
V of the velocity distribution n(v) must relate to a fixed frame. In the frame in which the 
Earth is at rest 



v± 



(14) 



where Ds is the distance to the source star, so, for fixed f^, v is a function of and to make 
further progress analytically, it is necessary to assume that both v and Vg have Gaussian 
distributions. Since we are working in the Earth's rest frame, these Gaussians will be centred 
on non-zero velocities v and v^, respectively. Consequently we have 



p(vs)ra(v) 



(27rasa) 



exp 



fv-v)2 



(15) 



2a2 2a2 

where ric is the volume density of clouds. Now 6 is the direction of v — {Dc/ Ds)vs, so with 
(p!^) we have 



v) 



v\ + 2t>_Lt>s COS 6 + vl, 



where 



-v« - V 



When we insert (^) into n(v) and integrate over 6 we find 



/ 



den(v) 



rir 



exp 



20-2 



(16) 



(17) 



(18) 



where /q is the usual modified Bessel function. Inserting this result in ( pTsD we find 



d^r 



dDAt Jo 



^R/t vl dv± 
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X 



d^Vs exp 



2^2 



(19) 



This result simplifies significantly if the sources have negligible random motions (i.e. asDc/Dg ^ 
a). In this case we can approximate the Gaussian distribution over Vg by 5^(vs — Vg), giving 



d^r 



^R/t dv^ vl 



^{2R/v^t) 
X exp 



2a2 



a' 



(20) 



with Vg replacing Vg in ( P7| ) . The rate of occultations integrated over all durations / dt d^V / dD^dt 
is readily obtained from either (|TU|) or (pUD: 



dr 



dDr -Kolo"^ 



X exp 
2nrR 



df I f ^ / d^v, 



(7^ 



2al 



dv± f _|_ exp 



2^2 



2a2 



V ^2 



(21) 



a" 



Since the central velocities Vg and v are relative to the Earth, equations (|T^ to (j2T|) yield 
occultation rates that depend on observing season. 

From equation (^Tj) T oc R/M (because Uc oc p/M), so for virialised clouds equation (^ 
indicates that the transit rate is controlled simply by the cloud temperature T, i.e. F oc T~^, 
independent of the scale size of the cloud. However, for fixed dimensionless impact parameter 
b/R the transit time t oc R from equation ([TT|) so dF/dt oc R^^T^^ when integrated over all 
impact parameters < < 1. 

In some applications the sources are widely distributed in distance D^. In this case we 
follow Kiraga & Paczyhski (1994) in evaluating weighted averages over of the rates (0) to 
(pip, with the weight factor taken to be UgD'^^'^^. Here rig is the source volume number density 
and (3 parameterizes the effects of the luminosity function upon the distance distribution 
of sources. In our calculations we adopt f3 = —2 for Dg > 100 pc and (3 = —1 at smaller 
distances. 



4 MODEL PREDICTIONS 

We explore two locations for potential cloud populations, the Galactic disc and halo, and 
we compute the rate of occultations for lines of sight to three targets: the Large Magellanic 
Cloud (LMC) at Galactic coordinates {l,b) = (280°, -33°), the Small Magellanic Cloud 
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(SMC) at (/, h) = (303°, -44°) and Baade's Window towards the Galactic Centre (GC) at 
(/, 6) = (1°,— 4°). These are all regions which have been intensively monitored by several 
experiments looking for microlensing events. As discussed in Section 0, a fourth possibility 
is the Andromeda Galaxy (M31), though in that case the formulae of the previous Section 
must be modified to take account of the fact that the sources in M31 are mostly unresolved. 
We therefore confine our attention in the remainder of this paper to the GC, LMC and SMC 
only. The LMC and SMC source stars are assumed to lie at fixed distances of 50kpc and 
60 kpc, respectively, whilst background sources towards Baade's Window are assumed to be 
distributed throughout the Galactic disc and bar. 

Since we are not concerned with the precise details of the halo structure, it is sufficient for 
our purposes to model the distribution of halo clouds as a simple softened isothermal sphere 
with a core radius of 5 kpc and a local density of 0.01 Mqpc"'^. The halo velocity distribution 
is taken to be Gaussian with a = 156 km s^^. For both clouds and sources in the Galactic 
disc we model their density distribution as a sech-squared profile with a scale length of 
2.5 kpc and scale height of 190 pc. We assume the clouds contribute up to a third of the disc 
density and so have a local density of 0.03 Mqpc"^. Strong local kinematical and dynamical 
constraints rule out the existence of a more significant population ( |Creze et al. 1998]) . We 
adopt a velocity dispersion of a = 25 km s~^ and a rotation speed v^ot = 220 km s~^ for both 
disc clouds and sources, independent of radius. 

For calculations towards the LMC and SMC we neglect random stellar motions in 
the clouds themselves though we do include the bulk motion of the Clouds (u, v, w) = 
(60, —155, 146) km s~^ ( [Jones, Klemola fc Lin 199^) , where the {u,v,w) coordinate system 
is centred on the Sun and the velocity components point, respectively, towards the Galactic 
Centre, the direction of rotation, and the North Galactic Pole. Equations ( pO]) and (^) are 
therefore sufficient to compute the halo cloud timescale distribution and rate towards these 
targets. For bulge sources we adopt the bar model of [Bissantz et al. (1997| ) with exponential 
scale length Om = 3.82 kpc, power-law scale length oq = 100 pc, axis ratios ">] = ( = 0.3 
and bar angle 0o = 25°. We assume a bar pattern speed of 70 km s~^ kpc~^ and a velocity 
dispersion a = 100 km s~^. Whilst our disc model is simpler than the double-exponential 
profile of pissantz et al. (1997|) its scale length is the same and the chosen scale height 
corresponds to the weighted mean of the scale heights in their model. Our overall disc-bulge 
mass normalization is therefore in good agreement. 

The motion of the observer is modeled as the sum of three distinct components: the 
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Figure 1. Seasonally-averaged differential (top) and cumulative (bottom) cloud transit timescale distributions for halo clouds, 
(solid lines) and disc clouds (dashed lines). The distributions are shown for three lines of sight: LMC (black); SMC (dark grey); 
GC (light grey). The figures assume the clouds have a radius R = 7 AU and virial temperature T = 10 K. The scaling for other 
values is dT/dt oc R-^T'^ and t oc R. 



motion of the LSR {u,v,w) = (0, Wrot, 0) km s~^; the Sun's motion in the direction of the 
Solar Apex [u, v, w) = (10, 5.2, 7.2) km s~^ ( [Binney fc IVlerritield 1998D ; and the motion of the 
Earth around the Sun, where for simphcity we adopt a circular orbit with % = 30 km s~^. 



4.1 Rate and timescale distribution 

Fig. shows the predicted distributions of occultation time for both halo and disc clouds 
towards the LMC, SMC and GC. The distributions are seasonal averages of the timescale 
distribution. We have adopted a radius R = 7 AU and virial temperature T = 10 K for 
both cloud populations, corresponding to a cloud mass of M = 10"'^ Mq from equation @. 
The top panel of Fig. |l] shows the differential distribution, whilst the lower panel gives the 
cumulative rate. For these parameters we see that the halo cloud distributions (solid lines) 
are peaked at shorter durations than for disc clouds (dashed lines), with the peak occurring 
around 60 days for all three target directions. For disc clouds the peak occurs at around 100 
days in the direction of the GC or around 300 days towards the LMC and SMC. 

For fixed cloud temperature the most detectable signals comes from disc clouds observed 
towards the GC and halo clouds towards the SMC, both with total rates of ~ 3000 events/yr 
per million stars. Halo clouds towards the SMC have the largest rate for transit durations 
below 500 days. Halo clouds observed towards the LMC and GC have smaller, but compara- 
ble, rates. These large rates in comparison with microlensing is due to the size of the clouds. 
For 1O~^M0 clouds their virial radius is 1.5-2 orders of magnitude larger than the typical 
Einstein radius of compact objects of the same mass, so the rate is correspondingly higher. 



10 Kerins, Binney & Silk 



01 



Gal. Centre 


Mag. Cloud 




~ ^ ^ observations 


observations y " 
/ 

/ 












/L^_^^ 












// 




\ ^ 


// 




LMC 

SMC 
GC 


// 





100 200 



Xg (degrees) 

Figure 2. Seasonal modulation of ttie occultation rate. Displayed as a function of Solar Ecliptic longitude, Aq, is the instanta- 
neous rate, F, expressed as a ratio of the seasonally-averaged rate, (F). Line coding is as for Fig. |^. The approximate observing 
season for the LMC and SMC is shown by the dark grey shaded region whilst the GC monitoring period is shown in light grey. 
The modulation curves are shown as thick lines during the relevant observing season and as thin lines otherwise. 

Both the timescales and the expected rate are easily within range of current microlensing 
experiments. Towards the LMC and SMC the disc cloud transit rate is lower by more than 
an order of magnitude at ~ 100 events/yr per million stars. Even when detection efficien- 
cies have been taken into account, microlensing data sets should contain hundreds or even 
thousands of transit events, so it is clear that much of the parameter space for both halo 
and disc cloud populations can be probed. 

4.2 Discriminating between disc and halo clouds 

Can we discriminate between halo and disc cloud populations? From Fig. |l| it is apparent 
that a rate of > 1000 events/yr per million sources towards the LMC and SMC would argue 
strongly for a halo origin. Another obvious diagnostic is the transit duration. Fig. |l| indicates 
that we should expect few halo cloud events with transit durations above 100 days, whilst 
more than three-quarters of the rate of disc clouds towards the GC comes from events with 
durations above 100 days. Of course larger halo clouds can produce longer transits so this 
is not a robust test. 

A potentially cleaner test is the seasonal modulation in the observed rate due the Earth's 
motion around the Sun. The modulation is most pronounced for nearby clouds, which will 
typically be disc clouds. Fig. ^ shows the modulation in the instantaneous transit rate F, 
as a ratio of the seasonally-averaged rate (F). The modulation is plotted as a function of 
Solar Ecliptic Longitude, Aq, where Aq = at the Vernal Equinox which occurs around 
March 21st. The dark grey region of the figure, spanning from September to March, indicates 
the six-month period when microlensing experiments monitor the Magellanic Clouds. For 
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the rest of the year the GC is observed, as shown by the hght grey region. Consequently, 
only half of the modulation period is accessible for any of the targets, so the modulation 
curves are shown as bold lines whilst the relevant target is observed and as thin lines when 
not. 

The largest modulation occurs for disc clouds observed towards the LMC and SMC. 
During the LMC/SMC observing season the rate swings from 0.8 to 1.2 times the seasonal 
average, a modulation of 20%. Assuming 100% detection efficiency, a 20% modulation could 
be detected at the 3 a level after four seasons of data for a million sources. Current ex- 
periments are monitoring closer to 10^ sources towards the LMC, though their detection 
efficiency to transits will be somewhat less than 100%. It's nonetheless clear that this mod- 
ulation should be detectable already. Halo clouds show only a small (few percent) variation 
in their rate towards the LMC and SMC, however their transit rate is much larger than for 
disc clouds. A sample upwards of 30000 transit events is required to detect a modulation of 
2% at greater than 3 a significance. Experiments monitoring 10'' sources would amass such 
a sample within four years, so even this small modulation should be evident within existing 
data sets. Towards the GC the modulation of the disc and halo cloud rate is similarly small 
but the rate is comparably large, so again the modulation should be apparent. In the case of 
halo clouds, observations towards the GC should show a maximum rate during mid-season, 
whilst disc clouds should exhibit a minimum. The modulation signal should therefore provide 
a good test for the origin of detected events. 

5 DISCUSSION 

Cold opaque clouds can be detected by searching data sets already compiled by groups 
searching for gravitational microlensing events. Rather than the flux excess that is the 
signature of microlensing, cold cloud events are characterized by a transient dimming of 
background sources. Like microlensing, this signal should be non-periodic and one should 
expect events to involve a representative population of background source stars. In selecting 
microlensing events one problem is to reject variable stars which may mimic a microlensing 
signal. The equivalent "background" for cloud transit events is potentially much smaller 
since there are fewer classes of objects which involve a signiflcant dimming of flux. One such 
class, eclipsing binaries, is unlikely to present much of a problem. For short-period systems 
their true nature would be evident over the lifetime of the surveys. Binary systems with 
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longer periods for which only one flux dropout is detected, should statistically comprise 
more massive stars, so we should expect the source stars of these "events" to be unrepre- 
sentative of the target stellar population as a whole, contrary to expectation for true cloud 
transit events. It is also likely that the light-curve of cloud transit events would have a form 
which generally would be inconsistent with that of eclipsing binaries. 

The previous section demonstrates that the transit rate and timescales of cold clouds 
in the disc and halo is well within the range of detectability if they constitute a significant 
population and are of a Jupiter mass scale. Whilst the total rate depends only on the 
temperature of the clouds, virialized clouds more massive than ~ 0.05 Mq would have 
transit durations exceeding the present baseline of microlensing experiments. 

The seasonal modulation of the rate provides a promising method to distinguish whether 
the clouds are in the disc or the halo. There are already strong arguments from the con- 



sideration of sub-mm sources ( [Lawrence 200 1|) that opaque cold clouds do not contribute 
much to the halo dark matter budget, though these arguments are sensitive to the precise 
dust properties of the clouds. Existing microlensing data sets represent a significant corpus 
of data which can provide an independent line of approach. Though originally conceived as 
a probe only of compact objects, microlensing surveys may also turn out to be one the most 
sensitive probes of diffuse cold clouds. 
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